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Abstract. The statistical properties of the return intervals between successive 1 -min volatilities of 30 liquid Chinese 
stocks exceeding a certain threshold q are carefully studied. The Kolmogorov-Smirnov (KS) test shows that 12 stocks 
exhibit scaling behaviors in the distributions of Tq for different thresholds q. Furthermore, the KS test and weighted 
KS test shows that the scaled return interval distributions of 6 stocks (out of the 12 stocks) can be nicely fitted by a 
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stretched exponential function /(t/t) 



-cv(t/t)'* 



with 7 ~ 0.31 under the significance level of 5%, where r is the 



mean return interval. The investigation of the conditional probability distribution Pq(r|ro) and the mean conditional 
return interval {t\to) demonstrates the existence of short-term correlation between successive return interval intervals. 
We further study the mean return interval {t\to) after a cluster of n intervals and the fluctuation F{1) using detrended 
fluctuation analysis and find that long-term memory also exists in the volatility return intervals. 

PACS. 89.65.Gh Economics; econophysics, financial markets, business and management - 89. 75. Da Systems obeying 
scaling laws - 05.45. Tp Time series analysis 



. 1 Introduction 

■ In recent years many concepts and methods from statistical 
' physics have been applied to the study of financial markets (l). 
. The statistical analysis of the waiting time between two suc- 
' cessive events has drawn much attention. Different definitions 
, of event refer to a variety of variables characterizing the prop- 
' erties of stock markets from certain aspects. For example, the 
[ persistence probability concerns the duration time that prices 
(or volatilities) always keep above or below their initial val- 
ues ilUSlIll, the exit time raises a concept similar to the per- 
sistence probability and focuses on the interval time between 
prices with fluctuations larger than a certain threshold [5 6,7^ 
[8ll9l [T0l . and the intertrade duration stands for the time interval 
between two successive independent trades 11 llll2l[T3l[T4l . 

Recently, extensive work has been done in the analysis of 
waiting time between extreme events in natural records, such as 
floods, temperatures and earthquakes Ill5lll6lfr7lll8l . The un- 
derstanding of the statistical properties of these extreme events 
is of great importance for the risk assessment of rare events. By 
investigating the return intervals between successive extreme 
events exceeding a certain threshold q, scaling behaviors are 
revealed in the return interval distributions for numerous com- 
plex systems. The dependence of the return interval distribution 
on the correlation of the original records has been carefully in- 
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vestigated Ill9ll20il2ni22ll23l . It shows that the scaling behav- 
ior of the return intervals might be due to the long-term mem- 
ory of the original records. It is well known that long memory 
is an important statistical property basically observed in stock 
volatilities yj. Therefore, one can expect that the scaling be- 
havior might also appear in the volatility records of stocks. 

A first effort was conducted by Yamasaki et ai, who used 
the daily data of US stocks to study the properties of the volatil- 
ity return intervals ll24l . Suppose that Y{t) is the price at time 
t, and the volatility is simply defined as the magnitude of loga- 
rithmic return R{t,St) = |hi(y(t)) - \n{Y{t - 6t))\. They 
studied the return intervals r between successive volatilities 
above a certain threshold q, and found that the distribution of 
volatility return intervals Pq{T) indeed obeys a scaling behav- 
ior. Wang et al. further studied the statistical properties of the 
volatility return intervals using the intraday data of US stocks 
Il25ll26ll27l . They found that Pq(r) exhibited a similar scal- 
ing behavior and the distribution could be approximated by 
a stretched exponential form similar to that found in climate 
records ifTTl . All these studies show that the volatility return 
intervals have long memory. Similar scaling behavior and long- 
term memory are observed in the daily and 1-min volatility 
return intervals of thousands Japanese stocks 1281 . Qiu et al. 
studied the volatility return intervals using the high-frequency 
intraday data of four liquid stocks traded in the Chinese market, 
which is an important emerging market, and found the return 
interval distribution also follows scaling behavior f29\. 
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Contrary to the scaling behavior found previously, Wang et 
al. analyzed the trade-by-trade data of 500 stocks which com- 
pose the S&P 500 index and found that the cumulative distri- 
bution of return intervals had systematic deviations from scal- 
ing and showed multiscaling behaviors 1301 . The analysis was 
based on the deduction that if the probability distribution obeys 
scaling behavior its cumulative distribution should also obeys 
scaling behavior They further studied the m-th moment of the 
scaled return intervals, and found the moment showed a certain 
trend with the mean interval which supports the finding that 
the return interval distribution exhibits a multiscaling behavior. 
This remarkable multiscaling behavior calls for a more careful 
scrutiny of the volatility return intervals by adopting solid sta- 
tistical methods rather than eyeballing. In addition, Lee et al. 
investigated the return intervals of 1-min volatility data of the 
Korean KOSPI index ISDl . They found that the interval return 
distribution had a power-law tail and no scaling was observed. 
However, it seems that they did not remove the intraday pattern 
from the intraday volatility series, which weakens their conclu- 
sion. Note that intraday pattern should be removed when deal- 
ing with intraday data, which will have heavy influence on the 
estimated return intervals for large thresholds. 

In this work, we use a nice high-frequency database |f32l 
to study the statistical properties of the interval returns. The 
two samples Kolmogorov-Smirnov (KS) test, which is a stan- 
dard way to check whether these two samples are drawn from a 
same distribution by comparing their cumulative distributions, 
is used to examine the coincidence of the interval distributions 
for different values of threshold q. Specifically, we study in de- 
tail 30 most liquid stocks to gain better statistics. Our results 
show that only 12 individual stocks pass the KS test and show 
scaling behaviors, and the remaining 18 stocks show nonscal- 
ing behaviors which seems consistent with that found by Wang 
et al. 1301 . Furthermore, a KS goodness of fit test is performed 
to study the particular form of the scaling function. The mem- 
ory effect is also investigated in this paper 

The paper is organized as follows. In Section|2l we explain 
the database analyzed and how the volatility return intervals 
are calculated. Section[3]studies the scaling and nonscaling be- 
haviors of the volatility return interval distributions using the 
KS tests. In Section H] we further study the memory effect of 
volatility return intervals. Section|5]concludes. 



2 Preprosessing the data sets 

Our analysis is based on the high-frequency intraday data of 
30 most liquid stocks traded on the Shanghai Stock Exchange 
and the Shenzhen Stock Exchange. These 30 stocks are most 
actively traded stocks representative in a variety of industry 
sectors, and thus have the largest sizes among all the stocks. 
The basic information about these stocks including the stock 
codes, company names and industrial sectors they belonging 
to is shown in Table [T] The dada record the trading prices ev- 
ery six to eight seconds from January 2004 to June 2006. The 
volatility is defined as the magnitude of logarithmic price return 
between two consecutive minutes, that is R{t) ~ \ \n{Y{t)) — 
hi(y(t — 1))|, where the price is the closest tick to a minute 
mark. Thus the sampling time is one minute, and the volatility 
data size is about 140,000 for each stock. 



For most stock markets, the intraday volatilities exhibit a 
U-shaped or L-shaped intraday pattern. This is so for the Chi- 
nese stocks l33l . When dealing with intraday data, this pat- 
tern should be removed I25,26,27l [30ll29l . Otherwise, the re- 
turn intervals distribution will exhibit daily periodicity for large 
thresholds. The intraday pattern A{s) is defined as 



A{s) 



N 



(1) 



which is the volatility at a specific moment s of the trading day 
averaged over all N trading days and i?'(s) is the volatility at 
time s of day i. The L-shaped intraday patterns of four typical 
stocks are illustrated in Figure [T] The intraday pattern exhibits 
a pronounced peak at the opening hours, and a minimum value 
just before the midday break of a trading day. In comparison 
with the western stock markets, the intraday pattern of the Chi- 
nese stock market surges soon after midday break which may 
due to the information aggregation during the midday break. 
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Fig. 1. (Color online) The 1-min interval intraday patterns of four rep- 
resentative stocks: 000002, 000625, 600036 and 600028. The curves 
are vertically shifted for clarity. 



To avoid the effect of this daily oscillation, the intraday pat- 
tern is removed by 

R{t) 



R'{t) = 



A{s) 



(2) 



Then we normalize the volatility by dividing its standard devi- 
ation 

r(i) = (3) 

We note that the raw 1-min Chinese stock returns follow Stu- 
dent distribution whose tail exponent is about 3.5 I34J . which 
ensures that the variance of R'{t) exists where the intraday pat- 
tern is removed. 

We study the return intervals t between successive volatili- 
ties 1 7' I exceeding a certain threshold q. The definition of return 
intervals is illustrated in Figure|2]for a typical stock. It is clear 
that the number of return intervals corresponding to a given 
threshold decreases with increasing threshold and the sum of 
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Table 1. Kolmogorov-Smimov test of the return interval distributions for q = 2 and 5 by comparing the statistic KS with the critical value 
CV at the 5% significance level. 



Code 


Stock name 


Industrial sector 




KS 




CV 


Scaling? 


000002 


China Vanke Co., Ltd 


Real estate 





0793 





0636 


No 


000503 


Searainbow Holding Co., Ltd 


Conglomerates 





0485 





0469 


No 


000625 


Chongqing Changan Automobile Co., Ltd 


Manufacturing 





0233 





0509 


Yes 


000839 


Citic Guoan Information Industry Co., Ltd 


IT 





0705 





0490 


No 


000858 


Wuliangye Yibin Co., Ltd 


Wine 





0803 





0500 


No 


000917 


Hunan TV&Broadcast Intermediary Co., Ltd 


Broadcasting 





0718 





0511 


No 


000930 


Anhui BBCA Biochemical Co., Ltd 


Biology 





0699 





0518 


No 


000983 


Shanxi Xishan Coal and Electricity Power Co., Ltd 


Power electronics 





0509 





0496 


No 


600000 


Shanghai Pudong Development Bank Co., Ltd 


Financials 





0699 





0538 


No 


600019 


Baoshan Iron & Steel CO., LTD 


Metals 





0223 





0672 


Yes 


600026 


China Shipping Development Co., Ltd 


Transportation 





0277 





0511 


Yes 


600028 


China Petroleum & Chemical Co. Ltd 


Petroleum & Chemistry 





0670 





0653 


No 


600030 


CITIC Securities Company Co., Ltd 


Financials 





0580 





0520 


No 


600036 


China Merchants Bank Co., Ltd 


Financials 





0436 





0576 


Yes 


600073 


Shanghai Maling Aquarius Co., Ltd 


Food 





0501 





0500 


No 


600088 


China Television Media Co., Ltd 


Broadcasting 





0540 





0497 


No 


600100 


Tsinghua Tongfang Co., Ltd 


IT 





0462 





0500 


Yes 


600104 


SAIC Motor CO., LTD 


Manufacturing 





0149 





0573 


Yes 


600110 


China-Kinwa High Technology Co., Ltd 


Manufacturing 





0911 





0506 


No 


600171 


Shanghai Belling Co., Ltd 


Electronics 





0377 





0530 


Yes 


600320 


Shanghai Zhenhua Port Machinery Co., Ltd 


Manufacturing 





0515 





0506 


No 


600428 


COSCO Shipping Co., Ltd 


Transportation 





0475 





0510 


Yes 


600550 


Baoding Tianwei Baobian Electric Co., Ltd 


Manufacturing 





0198 





0498 


Yes 


600601 


Founder Technology Group Co., Ltd 


IT 





0218 





0637 


Yes 


600602 


SVA Electron Co., Ltd 


Electronics 





0167 





0561 


Yes 


600688 


SINOPEC Shanghai Petrochemicl Co., Ltd 


Petroleum & Chemistry 





0421 





0566 


Yes 


600770 


Jiangsu Zongyi Co., Ltd 


Conglomerates 





0647 





0478 


No 


600797 


INSIGMA Technology Co., Ltd 


IT 





1179 





0647 


No 


600832 


Shanghai Oriental Pearl (Group) Co.Ltd 


Conglomerates 





0610 





0474 


No 


600900 


China Yangtze Power Co., Ltd 


Power electronics 





0574 





0556 


No 



3 Scaling and nonscaling behaviors of 
return interval distributions 

Several empirical studies show that the probability distribution 
function (PDF) of the return intervals obeys a scaling form: 

where f is the mean return interval which depends on the thresh- 
old q. The scaling form could be approximated by a stretched 
exponential function as 

f{x)^ce-"^\ (5) 

where c and a are two parameters and 7 is the correlation ex- 
ponent characterizing the long-term memory of volatilities. In 
this section, we investigate if these two properties hold for the 
30 Chinese stocks. 

3.1 Complementary cumulative distributions 

To make the observation of this possible scaling behavior (|4|l 
more clear, we study the complementary cumulative distribu- 
tion function (CCDF) of the scaled return intervals 

/•oo poo 
Cg{T/f) = / Pg{T)dT = / f{x)dx. (6) 

J T J t/t 



return intervals Tq are approximately identical to each other for 
different q. When the threshold q is very small, say less than 
the minimum of r{t), then all values equal to I. 




t [min] 

Fig. 2. Illustration of volatility return intervals for stock 000625, 
where r2, ra and r4 correspond to return intervals for g = 2, 3 and 4. 
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If the PDFs for different q obey the scaling form in Eq. (|4]i, the 
CCDFs of the scaled return intervals should also collapse onto 
a single curve. 

We calculate the CCDFs of the scaled return intervals of the 
chosen 30 stocks for a wide range of threshold q = 2, 3, 4, 5. 
Note that q is in units of standard deviations such that q = 5 
is a quite large volatility. The CCDFs of a representative stock 
000839 are shown in Figure |3j a). The curves for different val- 
ues of threshold q seem close to each other but do not col- 
lapse onto a single curve. A systematic deviation from scal- 
ing is clearly observed: the cumulative distribution decreases 
with the increase of threshold q for small scaled return inter- 
vals t/t < 1, while it increases with q for large scaled return 
intervals t/t > 1. This observation is similar to the multiscal- 
ing behavior reported by Wang et al. Il30l . 

However, we also find some stocks which apparently show 
scaling behaviors. In Figure |3jb), the CCDFs of a representa- 
tive stock 000625 that exhibits approximate scaling behavior 
are plotted. One observes that the curves for different values 
of threshold q approximately collapse onto a single curve. This 
indicates the return interval distributions of this stock obey a 
scaling behavior. 



3.2 Kolmogorov-Smirnov test of scaling in return 
interval distributions 

The eyeballing of complementary cumulative distributions of- 
fers a qualitative way of distinguishing scaling and nonscaling 
behaviors. Here we further adopt a quantitative approach, the 
Kolmogorov-Smirnov (KS) test. We use the KS test to compare 
two distributions for q = 2 and g = 5, which behave most dif- 
ferently among all the q value distributions. If the distributions 
for q = 2 and q = 5 pass the test, we can conclude that all the 
return interval distributions for different q values collapse onto 
a single curve and thus obey a scaling law. 

The standard KS test is designed to test the hypothesis that 
the distribution of the empirical data is equal to a particular dis- 
tribution by comparing their cumulative distribution functions 



(CDFs). Our hypothesis is that the two return interval distribu- 
tions for q ~ 2 and g = 5 do not differ at least in the com- 
mon region of the scaled return intervals 1281 . Suppose that 
F2 = 1 — C2 is the CDF of return intervals for q = 2 and 
F5 = 1 — C5 is the CDF of return intervals for q = 5. We 
calculate the KS statistic by comparing the two CDFs in the 
overlapping region: 

KS = ma^{\F2-F5\) . (7) 

When the KS statistic is smaller than a critical value denoted 
by CV (i.e., KS < CV), the hypothesis is accepted and we 
can assume that the distribution for g = 2 is coincident with 
the distribution for q = 5. The critical value at the significance 
level of 5% is CV = 1.36/ yjmn/ (m + n), where m and n are 
the numbers of interval samples for q = 2 and 5, respectively. 
In Table [T] is depicted the KS statistics and the corresponding 
critical values for the 30 stocks. One can see that 12 stocks 
including the stock shown in Figure |3lb) pass the KS test and 
consequently their return interval distributions obey scaling be- 
haviors, while other 18 stocks including the one shown in Fig- 
ure [3 a) do not pass the test. We find that stock 600030, which 
is the only stock also studied by Qiu et al. 1291 . does not pass 
the KS test at the significance level of 5%, but the statistic KS 
is very close to the critical value CV, i.e., 0.0580 ~ 0.0520. 
Further study shows that the stock 600030 passes the KS test 
at a significance level of 1.97% which indicates this stock may 
follow a relatively week scaling behavior. 

3.3 Kolmogorov-Smirnov test of the scaling function 

We have demonstrated that the return distributions of 12 stocks 
show scaling behaviors. To further study the particular form of 
the scaling function, we perform the KS goodness-of-fit test 
II35II36I . Empirical studies have shown that the scaling form 
could be approximated by a stretched exponential function as 
in Eq. (|5]l. In this case, our hypothesis is that the empirical 
distribution is coincident with the fitted stretched exponential 
distribution. Similar to the KS test we have conducted for two 
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empirical samples, we use the KS statistics to test whether the 
distributions for q = 2 and q = 5, which behave most differ- 
ently among all the q value distributions, are identical to a same 
fitted distribution for both q values in the overlapping region of 
the scaled return intervals. If both return interval distributions 
for (7 = 2 and g = 5 pass the test, we can conclude that all 
the return interval distributions for different q values follow a 
scaling function with stretched exponential form. 

The method is described as follows. Let Fq ~ 1 — Cq be the 
cumulative distribution of Tg, Fse the cumulative distribution 
from the fitted stretched exponential. The KS statistic defined 
in Eq. ^ becomes 



KS = max(|Fq - 



SE| 



9 = 2,5 



(8) 



We also use an variant of the KS statistic known as the weighted 
KS statistic, which is defined as follows |[36l 



KSW ^ max 



\/ -Fse(1 - ^se) 



(9) 



which is more sensitive on the edges of the cumulative distri- 
bution. In this case, the bootstrapping approach is adopted ll35l 
[36| . To do this, we first generate 1000 synthetic samples from 
the best fitted distribution and then reconstruct the cumulative 
distribution Fsi,„ of each simulated sample and its CDF Fsim.SE 
from integrating the fitted stretched exponential. We calculate 
the values of KS and KSW between the fitted CDF and the 
simulated CDF using 



KSsi 



F 



sim.SE 



I) 



and 



KSWsi 



\F^\y^ Fa 



sim,SE| 



\/^sim.SE(l — Fsim,SE) 



(10) 



(11) 



The p- value is determined by the frequency that KSsim > KS 
or KSWsiin > KSW. The tests are carried out for the 12 
stocks exhibiting scaling behaviors. The resultant p-values are 
depicted in Table |2] 

Consider the significance level of 1%. If at least one p- 
value for q = 2 or g = 5 of an individual stock is less than 
1%, then the null hypothesis that the empirical PDFs of this 
stock can be well fitted by a stretched exponential is rejected. 
According to Table |2] the null hypothesis is rejected for three 
stocks (600100, 600428, 600601) using the KS test and for two 
stocks (600428, 600601) using the KSW test. Under the signif- 
icant level of 5%, the null hypothesis is rejected for six stocks 
(000625, 600026, 600100, 600428, 600601, 600688) using the 
KS test and for five stocks (000625, 600026, 600100, 600428, 
600601) using the KSW test. We find that the KS test and the 
KSW test provide very similar results except that the KS test 
is slightly stronger ll36l . It is noteworthy to point out that the 
p-values of stocks 600602, 600019 and 600104 are very large, 
implying high goodness-of-fit of the stretched exponential to 
the empirical PDFs. In principle, the p-values of a stock are 
larger when the scaling of PDFs for different q is more sig- 
nificant, which is manifested by the fact that the statistic KS 
exceeds the critical value CV more in Table [T] 



Table 2. KS and KSW test of the return interval distributions for q — 
2 and 5 by comparing empirical data with the best fitted distribution 
and synthetic data with the best fitted distribution. 



Code 


9 


PR'S 


PKSW 


Code 


g 


Pks 


PKSW 


000625 


2 


0.224 


0.769 


600019 


2 


0.938 


0.742 




5 


0.045 


0.023 




5 


0.567 


0.574 


600026 


2 


0.010 


0.215 


600036 


2 


0.311 


0.432 




5 


0.021 


0.011 




5 


0.086 


0.098 


600100 


2 


0.003 


0.022 


600104 


2 


0.851 


0.973 




5 


0.157 


0.184 




5 


0.709 


0.866 


600171 


2 


0.158 


0.286 


600428 


2 


0.004 


0.001 




5 


0.401 


0.346 




5 


0.002 


0.007 


600550 


2 


0.226 


0.354 


600601 


2 










5 


0.059 


0.103 




5 


0.883 


0.130 


600602 


2 


0.995 


0.943 


600688 


2 


0.038 


0.386 




5 


0.996 


0.592 




5 


0.188 


0.056 



In all the fitting, the correlation exponent 7 for different 
stocks lies in a range [0.22, 0.44] with an average value 0.31, 
which is consistent with that observed in ifJTl . To show how 
good the stretched exponential fits the data, we illustrate in 
Figure |4] the empirical PDFs of Tq and the fitted curve for a 
representative stock 000625, which passes the KS and KSW 
goodness-of-fit tests at the significance level of 1 % but falls at 
the significance level of 5%. It is obviously that the curves for 
different values of q almost collapse onto a single curve which 
could be nicely fitted by a stretched exponential function. 
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Fig. 4. (Color online) Scaled probability distribution of scaled return 
intervals for different threshold q — 2, 3, 4, 5 for the representative 
stock 000625, which show scaling behavior. The solid curve is the 
fitted function ce"""" . 



4 Memory effect of return intervals 

The probability distribution may not fully characterize the prop- 
erties of volatility return intervals. The temporal correlation 
is known as another important observable independent of the 
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probability distribution. Empirical studies have revealed that 
the stock market volatilities are long-term correlated. We sup- 
pose that the long-term coiTelated volatilities would also af- 
fect the memory of the return intervals. Indeed, the memory 
effect is observed in the return intervals of US stock market 
and Japanese stock market. 



those shuffled return intervals, {t\to) fluctuates around a hor- 
izontal line close to 1 indicating that r is independent of the 
previous tq, as expected. This suggests that the memory of the 
return intervals may arise from the long-term memory of the 
volatilities. 



4.1 Short-term memory of return intervals 

To investigate the memory effect of the return intervals in Chi- 
nese stock market, we first calculate the conditional probabil- 
ity distribution Pg(T|ro), which is the probability to find an 
interval r immediately after the interval tq. Specifically, we 
study the conditional PDF for a bin of tq. The entire interval 
sequences are arranged in ascending order and partitioned to 
four bins with equal size. We illustrate the results using a typ- 
ical stock 000002. The results for all other 29 stocks are qual- 
itatively the same. Figure |5]plots Pg(r|ro) as a function of the 
scaled return intervals r/r for tq in the smallest and biggest 
bins. For large r/f , Pg(T|To) with tq in the biggest subset is 
larger than Pq{T\To) with tq in the smallest subset, and for 
smafl r/f, Pq(r|ro) with tq in the biggest subset is smaller 
than Pq{T\To) with tq in the smallest subset. This indicates 
that small intervals tq tend to be followed by small intervals 
T, while large intervals tq tend to be followed by large inter- 
vals T. Therefore, there exists short-term memory in the return 
intervals. 
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Fig. 6. (Color online) Scaled mean conditional return interval 
(r|ro)/f of scaled return intervals (filled symbols) and shuffled data 
(open symbols) for a typical stock 000002. The solid lines are guide- 
lines. 




Fig. 5. (Color online) Scaled conditional PDF Pq{T\To)f of scaled re- 
turn intervals r/f with ro in the largest 1/4 subset (open symbols) and 
smallest 1/4 subset (filled symbols) for a representative stock 000002. 
The solid lines are fitted curves. 



Calculating the mean conditional return interval (r|ro) im- 
mediately after tq is another alternative way to investigate the 
memory of return intervals. In fact, (r|ro) is the first moment of 
-fg('''|ro). In Figure|6]is shown (r|ro) of stock 000002, which 
exhibits a monotonously increasing tendency with the increase 
of the scaled return interval r/f. This also indicates that small 
r tends to follow small ro and big r tends to follow big ro. 
This result is consistent with that observed in Pg(r|ro). For 



4.2 Long-term memory of return intervals 




1234567 

cluster size n 



Fig. 7. (Color online) Scaled mean return interval (r|ro)/r after "+" 
cluster (filled symbols) and "-" cluster (open symbols) for a represen- 
tative stock 000002. 

To further study the long-term memory of return intervals, 
we investigate the mean return interval (r|ro) after a cluster of 
?? intervals that are all in a bin ro . The entire interval sequences 
are partitioned into two bins, separated by the median value 
of return intervals. Thus we investigate the mean return inter- 
val after a cluster of n intervals consecutively above and below 
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Fig. 8. (a) Detrended fluctuation F{1) of return interval T2 for four representative stocks. The curves are vertically shifted for clarity, (b) 
Extracted Hurst exponents H for large scales (circles) and small scales (squares) for the 30 stocks. 



the median value, denoted by "+" and "-" clusters. Figure |7] 
illustrates {t\tq) with respect to the cluster size n for a repre- 
sentative stock 000002. The results for other stocks are similar. 
In general, (t|to) after "+" cluster increases, while {t\to) after 
"-" cluster decreases as the cluster size n increases. In other 
words, a large mean return interval is likely followed by large 
consecutive return intervals and a small return interval is likely 
followed by small consecutive return intervals. This observa- 
tion shows the dependence of the return interval on its previous 
n consecutive return intervals. 

The detrended fluctuation analysis (DFA) method is applied 
to directly investigate the long-term correlation of return inter- 
vals. After removing a linear trend, the DFA method computes 
the detrended fluctuation F{1) of the time series within a win- 
dow of / points and determines the exponent H from the scal- 
ing form F{1) ~ , where H is the Hurst exponent II38II39I . 
For H > 0.5 the time series is long-term correlated, and for 
H = 0.5 the time series is uncorrected. The detrended fluctu- 
ation function F{1) of four representative stocks are shown in 
Figure [HIa). A crossover behavior is observed: F{1) for small 
scales I < 50 obeys a power law with an relatively small ex- 
ponent and F{1) for large scales I > 50 obeys a power law 
with a relatively large exponent. This phenomenon is also ob- 
served for other stocks, which is in line with other stock mar- 
kets 125)1 . Figure[8tb) illustrates the estimated Hurst exponents 
H of the 30 stocks. In both regions for small scales and large 
scales, the Hurst exponents of all the individual stocks are ap- 
parently greater than 0.5, which indicates the existence of long- 
term memory in return intervals. For shuffled data, the Hurst 
exponent displays a value close to 0.5. This confirms that the 
long-term memory may arise from the long-term memory of 
volatility records. 



5 Summary and conclusions 

In summary, we have studied the distribution and memory ef- 
fect of volatility return intervals for 30 most actively traded 
stocks on the Shanghai and Shenzhen Stock Exchanges. The 
Kolmogorov-Smirnov tests are performed to examine the scal- 
ing behavior of the return interval distributions as well as the 



particular form of the scaling function. We find that 12 stocks 
exhibit scaling behaviors in the volatility return interval distri- 
butions. We further find that the scaling form of 6 stocks having 
the best scaling can be excellently approximated by a stretched 
exponential function f{x) ~ e^"^ , and the mean value of the 
correlation exponent 7 is estimated to be around 0.31. 

The memory effect is also observed in the volatility return 
intervals of the 30 stocks. We first calculate the conditional 
PDF Pq{T\To) and mean conditional return interval (t|to), and 
find that the interval r is dependent of its closest previous inter- 
val To, which indicates the short-term memory of return inter- 
vals. We further investigate {t\tq) after a cluster of n intervals 
consecutively above and below the median value, and find that 
memory exists in a cluster of long return intervals. In addition, 
detrended fluctuation analysis of the return intervals shows that 
the return intervals have long-term memory. 

In a nutshell, our work focuses on the study of the statis- 
tical properties of volatility return intervals of Chinese stocks. 
Our main finding is that some stocks show scaling behaviors 
while the others do not. This novel feature shows the complex- 
ity of the Chinese stock market and help us better understand 
the properties of volatility return intervals. 
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gram for New Century Excellent Talents in University (No. NCET- 
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